Chaotic background of large-scale climate oscillations 

A. Bershadskii 
ICAR - P.O. Box 31155, Jerusalem 91000, Israel 



Abstract 

It is shown that the periodic alteration of night and day provides a chaotic 
dissipation mechanism for the North Atlantic (NAO) and Southern (SOI) cli- 
mate oscillations. The wavelet regression detrended daily NAO index for last 
60 years and daily SOI for last 20 years as well as an analytical continuation 
in the complex time domain were used for this purpose. 
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"And in the alteration of night and day, 
and the food We send down from the sky, 
and revives therewith the earth after its 
death, and in the turning about of winds 
are signs for people of understanding" . 

QURAN (45:4-5). 

1. Introduction 

It is now well known that climate is a nonlinear system. The nonlinear 
systems can exhibit a chaotic behavior. Both stochastic and deterministic 
processes can result in the broad-band part of the spectrum, but the decay 
in the spectral power is different for the two cases. The exponential decay 
indicates that the broad-band spectrum for these data arises from a deter- 
ministic rather than a stochastic process (cf. Figs. 1-4). For a wide class 
of deterministic systems a broad-band spectrum with exponential decay is a 
generic feature of their chaotic solutions (Ohtomo et al. 1995; Farmer 1982; 
Sigeti 1995; Frisch and Morf 1981). Also response of the chaotic systems to 
the periodic forcing does not always have the result that one might expect. 
Unlike linear systems, where periodic forcing leads to a periodic (peak-like) 
response, the chaotic response to periodic forcing can result in modification of 
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the rate of the exponential decay of the broad-band power spectrum. Thus, 
the chaotic response to periodic forcing can utihze the generic dissipation 
mechanism for the chaotic dissipative systems, especially in the case when 
the system's fundamental (internal) frequency is considerably lower than the 
frequency of the periodic forcing under consideration. 

The climate, where the chaotic behavior was discovered for the first time, 
is still one of the most challenging areas for the theory of chaos. The weather 
(time scales up to several weeks) chaotic behavior usually can be directly 
related to chaotic convection, while appearance of the chaotic properties for 
more long-term climate events is a non-trivial and challenging phenomenon. 

The solar day is a period of time during which the earth makes one rev- 
olution on its axis relative to the sun. During a part of the day the sun's 
direct rays are blocked (locally) by the earth. Therefore, the periodic daily 
variability of solar impact plays crucial role in high frequency climate behav- 
ior. It will be shown in present paper that just unusual properties of chaotic 
response to the daily periodicity of solar forcing provide an effective mech- 
anism for high frequency dissipation of the large-scale climate oscillations, 
presumably generated by the instabilities related to large scale topography. 

2. Large-scale climate oscillations 

One of the most significant and recurrent patterns of atmospheric vari- 
ability over the middle and high latitudes of the Northern Hemisphere is 
known as NAO - North Atlantic Oscillation. Climate variability from the 
subtropical Atlantic to the Arctic and from Siberia to the eastern boards 
of the North America is strongly related to the NAO (see, for a compre- 
hensive review Hurrell et al. (2003)). In site NAO (2010a) a projection of 
the daily 500mb height anomalies over the Northern Hemisphere onto the 
loading pattern (sec also NAO, 2010b) of the NAO was used in order to con- 
struct the daily NAO for last 60 years. Due to the natural climatic trends 
the daily NAO index time series is not a statistically stationary data set. In 
order to solve this problem a wavelet regression detrending method (Ogden, 
1997) was used in present investigation for the daily NAO time series. We 
used a symmlct regression of the data. Most of the regression methods are 
linear in responses. At the nonlinear nonpar ametric wavelet regression one 
chooses a relatively small number of wavelet coefficients to represent the un- 
derlying regression function. A threshold method is used to keep or kill the 
wavelet coefficients. In this case, in particular, the Universal (VisuShrink) 
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Figure 1: Spectrum of the wavelet regression detrended NAO data series (the NAO index 
data were taken from the site NAO, 2010a). In the semi-logarithmic scales used in the 
figure the straight line indicates the exponential decay Eq. (1). 

thresholding rule with a soft thresholding function was used. At the wavelet 
regression the demands to smoothness of the function being estimated are 
relaxed considerably in comparison to the traditional methods. Figure 1 
shows a spectrum of the wavelet regression detrended data (NAO, 2010a) 
calculated using the maximum entropy method. This method provides an 
optimal spectral resolution even for small data sets. The spectrum exhibits 
a broad-band behavior with exponential decay: 



A semi-logarithmic plot was used in Fig. 1 in order to show the exponential 
decay (at this plot the exponential decay corresponds to a straight line, for 
an explanation see below). 

Another example of the large-scale climate oscillations with the chaotic 
dissipation mechanism can be recognized in Tropical Atmosphere, where one 
of the most significant and recurrent patterns of atmospheric variability be- 
tween the western and eastern tropical Pacific is known as Southern Oscil- 
lation. The Southern Oscillation is related to the variabihty of the Walker 
circulation system: a circulation pattern characterized by sinking air above 
the eastern Pacific and rising air above the western Pacific. The Southern 
Oscillation is often considered as the atmospheric component of El Nino 
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Figure 2: Spectrum of the wavelet regression detrended SOI data series (the SOI data 
were taken from the site SOI, 2010). In the semi-logarithmic scales used in the figure the 
straight line indicates the exponential decay Eq. (1). 

phenomenon. The daily Southern Oscillation Index (SOI) has been calcu- 
lated based on the differences in air pressure anomaly between Tahiti and 
Darwin, Austraha (Troup, 1965) for last 20 years (SOI, 2010). Figure 2 
shows a spectrum of the wavelet regression detrended data calculated using 
the maximum entropy method. A semi-logarithmic plot was used in Fig. 2 
in order to show the exponential decay Eq. (1) (cf Fig.l). 

3. Chaos and exponential spectrum 

For a wide class of deterministic systems a broad-band spectrum with 
exponential decay is a generic feature of their chaotic solutions (Ohtomo et 
al. 1995; Farmer 1982; Sigeti 1995; Frisch and Morf 1981). Let us consider a 
relevant example. The Lorenz equations are given by: 

dx dy dz 

— — aiy — x), — = rx — y — xz, — — xy — bz (2) 
dt dt ^ ' dt ^ ^ ' 

The standard values producing a chaotic attractor are: a = 10.0, r = 
28.0, b = 8/3. Figure 3 shows a power spectrum for 2;-component of the 
Lorenz chaotic attractor generated by Eq. (2) (the spectrum was again cal- 
culated using the maximum entropy method). A semi- logarithmic plot was 
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Figure 3: As in Fig.l but for ^-component of the Lorenz chaotic attractor Eq. (2). In the 
semi-logarithmic scales used in the figure the straight line indicates the exponential decay. 

used again in Fig. 2 in order to show exponential decay (at this plot the 
exponential decay corresponds to a straight line). 

Another relevant example is a modulated Duffing oscillator. The driven 
Duffing oscillator has become a classic model for analysis of nonhnear phe- 
nomena and can exhibit both deterministic and chaotic behavior: Ott, 2002; 
Permann and Hamilton, 1992). It is well known that for strong nonlinearity 
a driven Duffing oscillator can exhibit chaotic properties and, in particular, 
exponentially decaying broad-band spectrum (see, for instance, Ohtomo et 
al. (1995)). In the case of a weak nonlinearity, however, one needs in an 
additional slow modulation of the driving force in order to obtain a chaotic 
behavior. In this separation between fast and slow motion is possible 

and in certain range of parameters the slow component exhibits chaotic prop- 
erties with an exponentially decaying broad-band spectrum (see, for instance. 
Miles (1984)). The driven Duffing oscillator is described by the equation 

X -\- ujqx(1 -\- Tjx^) + 2SuJox — A{t)cos{uJct) (3) 

where x denotes the temporal derivative of x, 5 is a damping parameter, r) 
is a nonlinearity parameter, cuq is the natural frequency for free oscillations, 
Uc is the carrier frequency. The amplitude A(t) can be slow modulated. 
Namely, A{t) = eujQCOsiuJrat + ^o); where ujm is a modulation frequency, 
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Figure 4: As in Fig. 3 but for slow component of the modulated Duffing oscillator. In the 
semi-logarithmic scales used in the figure the straight line indicates the exponential decay. 
The data were taken from Miles (1984). 



is a phase constant, e is a small constant amplitude. Except of a resonant 
domain: (w^ —Uq)/uJq and Wm/wo = 0{e'^^^) the nonlinearity is unimportant 
in the limit e — > 0. Following to Miles (1984) one can introduce a slow 
time T = ^auj. and pose the solution of Eq. (3) in the form (Van der Pol 
transformation) : 

x(t) = —[ui{T)cOs{ljJct) + U2{T)sin{ljJct)] (4) 
cr 

where Mi(t) and U2{t) are slowly varying amplitudes, a ~ 0(5^/^) is an 
arbitrary dimensionless positive constant. 

Then the equations for the slow components are 

^ = f/l(ni,M2), ^ = U2iUi,U2) + cosO (5) 

ar (IT 

where 9 = Qt + Oq, Q = 2uJm/o'UJc, and Ui{ui,U2) and f/2(wi,M2) are certain 
linear functions on ui and U2 (see for more details Miles (1984)). Numerical 
analysis of the equations (5) performed by Miles (1984) showed that for cer- 
tain values of the parameters the slow component of the oscillations exhibits 
chaotic properties. In particular, the broad-band spectra with exponential 
decay were observed by Miles (1984) for these values of parameters. Fig- 
ure 4 shows an example of such spectrum (f2 = 0.75, the data were taken 
from Miles (1984)). The spectra were calculated by Miles (1984) using a 
fast- Fourier-transform with frequency / i — > r. One can see a clear peak 
corresponding to the fundamental frequency = VL/2tt. We have drawn 
the straight line in Fig. 4 in order to indicate the exponential decay: Eqs. 
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(1). The same exponential decay - Eq. (1), can be observed for all 'chaotic' 
values of the parameter Q, i.e. it has an universal nature. 

Nature of the exponential decay of the power spectra of the chaotic sys- 
tems is still an unsolved mathematical problem. A progress in solution of this 
problem has been achieved by the use of the analytical continuation of the 
equations in the complex domain (see, for instance, Prisch and Morf (1981)). 
In this approach the exponential decay of chaotic spectrum is related to a 
singularity in the plane of complex time, which lies nearest to the real axis. 
Distance between this singularity and the real axis determines the rate of the 
exponential decay. For many interesting cases chaotic solutions are analytic 
in a finite strip around the real time axis. This takes place, for instance for 
attractors bounded in the real domain (the Lorenz attractor, for instance). 
In this case the radius of convergence of the Taylor series is also bounded 
(uniformly) at any real time. Let us consider, for simplicity, solution u{t) 
with simple poles only, and to define the Fourier transform as follows 

u(uj) = (27r)-^/2 dr e-'^^'uiT) (6) 

J-Te/2 

where the function ■u(r) is defined on the interval — Te/2 < r < Te/2 with 
periodic boundary conditions. Then using the theorem of residues 

u{oj) — i{27ry^'^ ^ Rj eyi-p(iojXj — \ojyj\) (7) 
j 

where Rj are the poles residue and Xj + iyj are their location in the rele- 
vant half plane of the complex time, one obtains asymptotic behavior of the 
spectrum E{(x!) = \u{(x!)\'^ at large a; 

E{f) ~ exp(-47r/ y„,„) (8) 

where ou — 27r/ and Umin is the imaginary part of the location of the pole 
which lies nearest to the real axis. 

If in the considered case ymin = Id, then we obtain the exponential decay 
shown in Fig. 1 (cf Eq. (1)). In order to understand the reason for ymin = Id 
we need in a nonlinear dynamic climate model where the daily periodic solar 
impact results in such position of the pole nearest to the real axis. But 
even before constructing such dynamic model it is already clear that just the 
daily periodicity of the solar impact is responsible for the chaotic exponential 
decay shown in Fig. 1. 
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Figure 5: The pertaining average maximal Lyapunov exponent at the pertaining time, 
calculated for the same data as those used for calculation of the spectrum (Fig. 1). The 
dashed straight line indicates convergence to a positive value. 

A confusion can arise due to the 1 day period also being the samphng fre- 
quency for the NAO and SOI time series. However, the poles location in the 
complex time plane is presumably determined by the system's dynamics and 
not by the sampling rate. For the considered above examples, for instance, 
the exponential decay rate (Figs. 3 and 4) is not equal to Atts, s being the 
sampling rate for the corresponding time series, as it would be otherwise. 

Additionally to the exponential spectrum, let us check the chaotic charac- 
ter of the wavelet regression detrended NAO data set calculating the largest 
Lyapunov exponent : Xmax- A strong indicator for the presence of chaos in 
the examined time series is condition Xmax > 0. If this is the case, then we 
have so-called exponential instability. Namely, two arbitrary close trajecto- 
ries of the system will diverge apart exponentially, that is the hallmark of 
chaos. To calculate Xmax we used a direct algorithm developed by Wolf et al. 
(1985), Kodba et al. (2003). Figure 4 shows the pertaining average maximal 
Lyapunov exponent at the pertaining time, calculated for the same data as 
those used for calculation of the spectrum (Fig. 1). The largest Lyapunov 
exponent converges very well to a positive value Xmax — 0.053(i^^ > 0. 
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4. Conclusions 

The exponential decay of the broad-band spectrum of the wavelet regres- 
sion detrended NAO and SOI time series (Figs. 1 and 2) indicates a chaotic 
nature of the high-frequency components of these climate oscillations. In Fig. 
1 one can readily recognize a peak corresponding to a dominant frequency 
of internal variability (corresponding to a period approximately equal to 40 
days, cf. with the fundamental frequency of the chaotic systems: Figs. 3 
and 4). It should be noted, that the 40-day oscillations are well known as 
an intrinsic mode of the Northern Hemisphere extratropics, which is exited 
presumably due to instabilities related to large-scale topography (see, for in- 
stance, Magana 1993; Marcus, Ghil and Dickey 1994; Ghil and Robertson 
2002). In this case the dominant frequency of internal variability (corre- 
sponding to the period ~40d) is considerably smaller than the daily solar 
forcing frequency, which determines the rate of the spectral exponential de- 
cay (Eq. (8)). Thus, the exponential chaotic decay (Fig. 1) can be considered 
as a dissipation mechanism for these oscillations. Figure 2 shows analogous 
picture for the SOI data. 
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